Stirring trapped atoms into fractional quantum Hall puddles 
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We theoretically explore the generation of few-body analogs of fractional quantum Hall states. 
We consider an array of identical few-atom clusters (n = 2,3,4), each cluster trapped at the node 
of an optical lattice. By temporally varying the amplitude and phase of the trapping lasers, one 
can introduce a rotating deformation at each site. We analyze protocols for coherently transferring 
ground state clusters into highly correlated states, producing theoretical fidelities (probability of 
reaching the target state) in excess of 99%. 

PACS numbers: 67.85.-d,73.43.-f,37.90.+j,03.75.Lm 



Cold atom experiments promise to produce unique 
states of matter, allowing controllable exploration of ex- 
otic physics. For example, since rotation couples to neu- 
tral atoms in the same way that a uniform magnetic field 
couples to charged particles, many groups are excited 
about the possibility of producing analogs of fractional 
quantum hall states [3, 0j S 0| • In particular, if a two di- 
mensional harmonically trapped gas of bosons is rotated 
at a frequency Q sufficiently close to the trapping fre- 
quency a;, then the ground state will have vortices bound 
to the atoms - an analog of the binding of flux tubes 
to electrons in the fractional quantum hall effect. The 
ground state will be topologically ordered and possess 
fractional excitations. Technically, the difficulty with re- 
alizing this goal experimentally has been that it requires 
O to be tuned to a precision which scales as 1/n, where 
n is the number of particles. Responding to this imped- 
iment, several authors 0, H, H| have proposed studying 
clusters with n < 10. Such puddles possess many of the 
features of a bulk quantum hall liquid, and producing 
them would be a great achievement. Here we propose and 
study protocols for producing strongly correlated clusters 
of rotating atoms. 

The issue prompting this investigation is that in such 
clusters there are very few mechanisms for dissipating 
energy, and hence experimentally producing the ground 
state of a rotating cluster is nontrivial. First, the small 
number of particles results in a discrete spectrum, and 
leaves few kinetic paths. Second, in the strongly cor- 
related states of interest the atoms largely avoid each 
other, further blocking the kinetics. On these grounds, 
one should not expect to be able to cool into the ground 
state. Instead we advocate a dynamical process where 
one coherently drives the system into the strongly cor- 
related state through a well-planned sequence of rotat- 
ing trap deformations. This approach is based upon an 
analogy between the states of these atomic clusters, and 
the energy levels of a molecule. By deforming the har- 
monic trap, and rotating the deformation, one couples 
the many-body states in much the same way that an 
oscillating electric field from a laser couples molecular 
states. We consider a number of pulse sequences, find- 
ing that one can rapidly transfer atoms to a strongly 



correlated state with nearly unit efficiency Following a 
proposal by Popp et al. [3[ , experimentalists at Stanford 
have achieved considerable success with a related pro- 
cedure, where one slowly increases the rate of rotation, 
adiabatically transfering bosonic atoms from an initially 
non-rotating state, to an analog of the Laughlin state [f|. 
One could also imagine implementing more sophisticated 
protocols such as rapid adiabatic passage |6j. 

To achieve sufficient signal to noise, any experimental 
attempt to study small clusters of atoms must employ 
an ensemble of identical systems: for example by trap- 
ping small numbers of atoms at the nodes of an optical 
lattice. When formed by sufficiently intense lasers, this 
lattice will isolate the individual clusters, preventing any 
"hopping" from one node to another. We will not ad- 
dress the very interesting question of what would happen 
if the barriers separating the clusters were lowered. By 
using Altering techniques, one can ensure that the same 
number of atoms sit at each node ||. A rotating defor- 
mation of each microtrap can be engineered through a 
number of techniques. For example, if the intensity of 
the lattice beams forming a triangular lattice are mod- 
ulated in sequence, then a rotating quadrupolar defor- 
mation is be produced. A more versatile technique is to 
modulate the phases between counterpropagating lattice 
beams. Changing these phases uniformly translates the 
lattice sites. If one moves the lattice sites around faster 
than the characteristic times of atomic motion (10 _5 s) 
but slow compared to the times for electronic excitations 
(10~ 15 s) then the atoms see a time averaged potential. 
This technique, which is closely related to the time or- 
bital potential traps pioneered at JILA Q , can produce 
almost arbitrary time dependent deformations of the in- 
dividual traps which each of the clusters experiences [f| • 
Each cluster feels the same potential. 

Once created, the ensemble of clusters can be experi- 
mentally studied by a number of means. In situ probes 
such as photoassociation Q and RF spectroscopy 0] re- 
veal details about the interparticle correlations. In the 
regime of interest, time-of-flight expansion, followed by 
imaging, spatially resolves the ensemble averaged pre- 
expansion density. This result follows from the scaling 
form of the dynamics of lowest Landau level wavefunc- 
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tions |4J. 

We model a single cluster as a small number of two- 
dimensional harmonically trapped bosonic atoms. The 
two dimensionality can be ensured by increasing the in- 
tensity of the lattice beams in the perpendicular direc- 
tion. Neglecting the zero-point energy, one finds that in 
the frame rotating with frequency SI, the single parti- 
cle harmonic oscillator eigenstates have the form Ejk = 
h(uj — tt)k + K(u> + Sl)j; (j, k = 0, 1, . . .). In typical lat- 
tices, the interaction energy U/h ~10 kHz is small com- 
pared to the small oscillation frequency u ~100 kHz [Icj |. 
Therefore the many-body state will be made up of single 
particle states with j = 0: the lowest Landau level, with 
wavefunctions of the form 
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where z — (x + iy)/d with d = yJh/Mui is the complex 
representation of the coordinate in the plane measured 
in units of the oscillator length, where M is the atomic 
mass. Including interactions, the many-body Hamilto- 
nian for a single cluster is then 

Hlll = ^^iHuj - n)a]aj + ^ Vjki m a]alaia m (2) 

j jklm 

where a m is the annihilation operator for the single parti- 
cle state ip m . For point interactions the interaction kernel 
is 



Vjklm — ~2°3+k-l-m* 
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where U = \/2/7r h 2 a/(Md z d 2 ) is the on-site interaction 
between two particles in the same well, a is the three di- 
mensional s-wave scattering length and d z is the oscillator 
length in the transverse direction. As has been explored 
in depth by previous authors [l|, 0, H, H[ , for a given total 
number of particles n, and angular momentum projec- 
tion L, the Hamiltonian @ is a finite matrix which is 
readily diagonalized. Example spectra are shown in fig- 
ures Q] through [3l We plot the spectra as energy versus 
angular momentum, with Q = u. Spectra at other ro- 
tation speeds are readily found by "tilting" the graphs - 
the energy of a state with angular momentum projection 
L is simply shifted up by h(u> — Q)L. 

We imagine applying to each cluster a rotating single 
particle potential (in the lab frame) of the form 



H s (r,t) = V p (t) [z 
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where m determines the symmetry of the deformation 
(e.g. m = 2 is a quadrupolar deformation), the enve- 
lope function V p {t) is the time-dependent amplitude of 
the deformation, and Sl p is the frequency at which the 
perturbation rotates. We will mainly focus on the case 
m = 2. When restricted to the lowest Landau level, this 
potential generates a coupling between the many-body 
states which in the rotating frame is expressed as 
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FIG. 1: (Color Online) Transfering small clusters from non- 
rotating ground state to v = 1/2 Laughlin state using ro- 
tating quadrupolar (m = 2) deformations. Left: Interac- 
tion energy (in units of 17/2) of quantum states of harmon- 
ically trapped two dimensional clusters as a function of to- 
tal angular momentum projection L in units of h. Excita- 
tion paths are shown by arrows. Central: squared overlap 
(fidelity) of \ip(t)) with the initial (solid) and final (dashed) 
states as a function of the duration of a square pulse. Right: 
Fidelities as a function of time for an optimized Gaussian 
pulse of the form e - '* - * - 1 ^ T . Time is measured in units of 
to = h/U ~ 10 -4 s. For n = 2, the peak perturbation am- 
plitude is V p = 0.05(17/2), lj - Q p = 2.0(17/2), and a Gaus- 
sian pulse time of r = 24ro. For n — 3, r = 102ro and 
cu - tt p = 2.046(17/2) and 2.055(?7/2) for the Gaussian and 
square cases, respectively. For n = 3, nonlinear effects (cou- 
pling with near-resonant levels) shifted the optimal frequency 
away from the linear response expectation, uj — fi p = 2(U/2). 



with vi m = 2- m / 2 (l + m)\/y/l\(l + m)l As such it only 
couples states whose total angular momentum projection 
differs by m. For our calculation will work in the co- 
rotating frame with SI = f2 p , where the only time depen- 
dence is given by V p (t). 

We wish to implement a 7r-pulse, where the ampli- 
tude V p (t) is engineered so that after the pulse, a clus- 
ter is transfered from its initial state to a target state 
of our choosing. If the perturbation coupled only two 
states, this would be a straightforward procedure. The 
frequency fi p is selected so that the initial and target 
state are degenerate in the rotating frame. For any fi- 
nite Vp, the system Rabi flops between the two cou- 
pled states, and by turning off the perturbation at the 
right time one ends up in the target state with unit 
probability. The present example is more complicated, 
as there are many states coupled by the perturbation. 
The basic idea however remains sound: one still chooses 
Sip to make the initial and final state degenerate. The 
time dependence of V p (t) should be tailored to mini- 
mize the coupling to unwanted states. These stray cou- 
plings could be particularly disasterous, because the cou- 
pling between the initial and target state are generi- 
cally quite high order m V p . As a particularly rele- 
vant example, we consider transfering clusters from the 



3 



ground state (with L — 0) to the the v = 1/2 Laughlin 
state ip L (zi, ■■■z n ) = rii<j(^-^) 2ex P(- Ej MV^I), 
which has angular momentum n(n — 1). Using a per- 
turbation with m = 2, this requires a n(n — l)/2-order 
process. A picturesque way of thinking about the dy- 
namics in the presence of the perturbation is to map 
the problem onto the motion of a particle on a compli- 
cated "lattice". The states of the unperturbed system 
are analogous to "lattice sites", while the perturbation 
produces a "hopping" between sites. The goal is to en- 
gineer a time-dependent hopping which efficiently moves 
the "particle" from a known starting position to a desired 
ending position. The transfer efficiency is measured by 
the probability that the system is in the target state i/jt 
at the end of the time evolution: we plot this probability 
- known as the fidelity - as a function of time, given by 

/(t) = |(^ T |^)>| 2 - 

As this analogy emphasizes, the problem of trans- 
fering a quantum system from one state to another is 
generic. Muller, Chiow, and Chu [11| recently consid- 
ered how one can optimize pulse shapes to produce high 
order Bragg diffraction, while avoiding transfering atoms 
into unwanted momentum states. These authors devel- 
oped a formalism for calculating the fidelity by adiabat- 
ically eliminating the off-resonant states. They found 
that Gaussian pulse shapes greatly outperformed simple 
square pulses. This result is natural, as the smoother 
pulses have a much smaller bandwidth. 

We numerically solve the time dependent Schrddinger 
equation, truncating our Hilbert space at finite total an- 
gular momentum L — n(n — 1) + 4 for n = 2, 3 and 
L = n{n — 1) + 8 for n — 4. We have numerically verified 
that changing this cutoff to higher values has negligible 
effects. Figure [T] shows / in the case of n = 2 and n = 3 
for square and Gaussian pulses. For the square pulse the 
fidelity is shown as a function of pulse length. For the 
Gaussian pulse, a fixed pulse duration is used, and the 
fidelity is shown as a function of time. For n = 2, where 
only two states are involved, the pulse shape is irrele- 
vant. For n — 3, where there is a near-resonant state 
with L = 4, the Gaussian pulse shape greatly outper- 
forms the square pulse, producing nearly 100% transfer 
efficiency in 10's of ms, even for a very weak perturbation. 

For n > 3 we find that these high order processes be- 
come inefficient. For m = 2 the coupling between the ini- 
tial and final state scale as (V^/Z7) n ( n-1 )/ 2 , making trans- 
fer times unrealistical long unless one drives the system 
into a highly nonlinear regime. As illustrated in figure [2] 
this difficulty can be mitigated by using perturbations 
with higher m. There, for illustration, we consider excit- 
ing a 3-particle cluster from the lowest energy L — 2 state 
to the L = 6 Laughlin state. The second order m = 2 
pulse requires much longer than the first order m = 4 
pulse. An interesting aside is that one would naively 
have expected that the resonant I = 4 state would make 
the second-order process extremely inefficient. It turns 
out that the coupling to that state is fortuitously zero. 

Further improved scaling can be arranged by using a 



sequence of ir pulses. One transfers the cluster from one 
long-lived state to another. Since the number of pulses 
scales as the angular momentum, the transfer time is 
then quadratic in the angular momentum, rather than 
exponential. One can also tailor the path to maximize 
the fidelity of each step. Some two-pulse sequences are 
shown in figure [31 The guiding principle in designing the 
pulse sequences is that in each step one wants as few as 
possible near-resonant intermediate states. 
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FIG. 2: (Color Online) Using a rotating m-fold symmetric 
perturbation to drive n = 3 particle clusters from L — 2 
to the v = 1/2 Laughlin state. Left:path on the energy 
level diagram. Center: second-order process coming from 
a deformation with m — 2. Right: direct transition pro- 
duced with m = 4. Solid (dashed) lines are fidelities with 
the initial (Laughlin) state. In both cases the peak deforma- 
tion is V p = 0.05(17/2). Both use a Gaussian pulse. The 
frequencies and pulse times r we used for m — 2, 4 were 
u-tlp = (3.00/2) (U/2), 3.035(172) and t/t = 218, 21. Note 
how much more rapid the direct process is. 

Summary. — We have shown it is possible to use time 
dependent trap perturbations to coherently transfer bo- 
son clusters from nonrotating ground states to analogs 
of fractional quantum hall states. We achieve fidelity 
/ > 99% for n — 2,3 using very weak rotating m = 2 
deformations, whose duration is of order tens of ms. Us- 
ing a two-pulse sequence, we achieve similar results for 
n = 4. We find that smooth Gaussian pulses are much 
more effective than square pulses, and that further effi- 
ciency can be gained by using higher order perturbations 
of the form z m with m > 2. 

We briefly compare our technique with Ref. l 3|'s pro- 
posal. While our approaches share the use of a rotating 
time-averaged optical lattice potential, our proposal of- 
fers significant differences and advantages. While Ref. [|[ 
suggests an adiabatic evolution, we propose a coher- 
ent evolution - analogous to a Rabi oscillation - to the 
Laughlin state. This has the advantage of being faster, 
easier to implement, and more robust. For a slightly 
smaller perturbation relative to the adiabatic method, we 
achieve fidelity ~ 1 in contrast to the adiabatic method's 
0.97 fidelity. Moreover, our method requires half the 
time. More importantly, the adiabatic method requires 
carefully navigating a path through possible rotating po- 
tential strengths and frequencies as a function of time. In 
contrast, our method requires only setting the pulse dura- 
tion and strength, and is thus more easily implementable 
and less susceptible to small experimental errors. 

This technique will allow the efficient creation of 
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FIG. 3: (Color Online) Transfering atoms using multiple 
pulses. Left: paths from initial to Laughlin states for n = 3,4. 
Right: Solid line is the fidelity with the initial state, dot- 
ted with the intermediate (L,E) — (2h,3(U/2)) state, and 
the dashed line with the Laughlin state. All pulses are 
Gaussians. Despite using multiple pulses, this technique is 
faster than using a higher order m = 2 pulse. The frequen- 
cies shape (m), and pulse times (r) for the N = 3 
sequence were h(u - D. p )/{U/2) = 3.00,3.035, m = 2,4, 
and t/tq = 16.95, 19.2. For both, V p = 0.05(t//2). For 
N = 4, using two pulses with m = 2 and Vp = 0.2([//2), 
we achieve > 98% fidelity after a total two-pulse sequences 
with h(u - Q p )/(U/2) = 3.130, 1.0376 and r/r = 82.5, 87.0. 



bosonic quantum Hall puddles - a state of matter which 
has not yet been observed. The clusters produced will 
be orbitally entangled, have strong interparticle correla- 
tions, have fractional excitations, and possess topological 
orders 0, 0, Hi- Although current experiments, and 
the present theory, is focussed on the small- atom limit, 
it would be exciting to apply these techniques to larger 
collections of atoms, producing true analogs of fractional 
quantum Hall states. The main difficulty is that the spec- 
tra become dense as n increases, requiring one to set fl p 
to extremely high precision. By carefully choosing the 
trajectory, taking advantage of gaps in the spectrum, one 
might be able to overcome such difficulties. 



Finally, we mention that our approach allows one to 
drive the system into almost arbitrary excited states. 
This may, for example, be important for using quantum 
hall puddles in a topological quantum computing scheme 

EH- 
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